Direct evidence of a zigzag spin chain structure in the honeycomb lattice: A neutron 
and x-ray diffraction investigation on single crystal Na2lr03 
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We have combined single crystal neutron and x-ray diffractions to investigate the magnetic 
and crystal structures of the honeycomb lattice Na2lr03. The system orders magnetically below 
18.1(2) K with Ir 4+ ions forming zigzag spin chains within the layered honeycomb network with 
ordered moment of 0.22(1) fis/lr site. Such a configuration sharply contrasts the Neel or stripe 
states proposed in the Kitaev-Heisenberg model. The structure refinement reveals that the Ir atoms 
form nearly ideal 2D honeycomb lattice while the IrOe octahedra experience a trigonal distortion 
that is critical to the ground state. The results of this study provide much-needed experimental 
insights into the magnetic and crystal structure crucial to the understanding of the exotic magnetic 
order and possible topological characteristics in the 5d-electron based honeycomb lattice. 

PACS numbers: 61.05.cf,75.25.-j,71.70.Ej 



The 5d-based iridates have recently become a fertile 
yet largely uncharted ground for studies of new physics 
driven by the spin-orbit coupling (SOC). It is now rec- 
ognized that the SOC (0.4 - 1 eV), which is proportional 
to Z 4 (Z is the atomic number), plays a critical role in 
the iridates, and rigorously competes with other relevant 
energies, particularly the on-site Coulomb interaction U 
(0.4 - 2.5 eV), which is significantly reduced because of 
the extended nature of the 5d orbitals. A new balance be- 
tween the competing energies is therefore established in 
the iridates and drives exotic states seldom or never seen 
in other materials. Recent experimental observations and 
theoretical proposals for the iridates has already cap- 
tured the intriguing physics driven by SOC: J e g = 1/2 
Mott state jJ-Qj spin liquid in hyper-kagome structure 
0, high-Tc superconductivity Weyl semimetal with 
Fermi arcs IE | , correlated topological insulator with large 
gaps 0, GJ, Kitaev mode 0, 3D spin liquid with 
Fermionic spinons fl3j |. etc. 
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Of all iridates studied so far, Na2lr03 
has inspired a great deal of experimental and theoreti- 
cal efforts. In essence, the honeycomb lattice Na2lr03 is 
predicted be a topological insulator or a layered quan- 
tum spin Hall insulator 10]. However, conspicuous dis- 
crepancies among various theoretical proposals and ex- 
perimental observations clearly point out a lack of the 
much-needed characterization of the magnetic and crys- 
tal structures of Na 2 Ir03, whose band topology could 
vary significantly with slight variations in the crystal 
structure. This situation chiefly originates from the fact 
that the heavy transition metals such as Ir strongly ab- 
sorb neutrons, which makes a comprehensive neutron 
study on the single crystal a nontrivial challenge. 



In this letter, we report a combined study of neu- 
tron scattering and x-ray diffraction on relatively large, 
thin single-crystal Na2lr03. This study reveals that Ir 4+ 
ions order magnetically below 18.1(2) K, and that they 
form zigzag spin chains along the a axis of the honey- 
comb structure with an ordered moment of 0.22(1) /is/Ir. 
Moreover, the structural refinements illustrate that the Ir 
atoms feature a nearly perfect two-dimensional (2D) hon- 
eycomb lattice and a trigonal distortion characterized by 
the IrOg octahedra deviating from a high-symmetric cu- 
bic environment. These results are vastly different from 
those of the previously reported powder x-ray diffraction 
study, where the honeycomb was characterized by three 
distinct bond lengths [14| . 

Single crystals of Na2lr03 were grown using a self-flux 
method from off-stoichiomctric quantities of Ir02 and 
Na2C03. Technical details are described elsewhere 0- 
6]. The crystals have circular basal area corresponding 
to the honeycomb plane with diameters of ~ 10 mm and 
thickness ~ 0.1 mm [see inset of Fig. 1]. Such geome- 
try provides a unique advantage to significantly allevi- 
ate the technical difficulty due to the inherent neutron 
absorption of the iridates. Energy dispersive x-ray spec- 
troscopy (EDX) using a Hitachi /Oxford SEM/EDX indi- 
cates a perfect stoichiometry of Na2lr03 throughout the 
crystals studied. 

The x-ray diffraction measurements were performed 
using a Bruker SMART APEX CCD diffractometer with 
Mo K Q radiation. More than 40 crystals from four dif- 
ferent growth runs were screened at 125 K and full data 
sets were collected on four crystals. The neutron diffrac- 
tion measurements were carried out at the HB1A triple 
axis spectrometer and HB3A four circle diffractometer 



2 



at the High Flux Isotope Reactor at the Oak Ridge Na- 
tional Laboratory with fixed incident neutron wavelength 
of A = 2.367 A and 1.536 A, respectively. For the HB1A 
diffraction measurement, the crystal is aligned in the 
(0, k, I) scattering plane to allow the probing of various 
magnetic reflections. 






FIG. 1. (Color online) (a) Crystal structure of Na2lr03 with 
C2/m symmetry, (b) The honeycomb lattice formed by Ir 
atoms within the basal plane with nearly equivalent distance 
between neighboring Ir atoms. The dashed line denotes the 
unit cell, (c) Specific heat C(T) for Na2lr03 at H = and 
H = 14 T. Note that the application of a magnetic field of 
14 Tesla suppresses the transition temperature by only 0.5 K, 
apparently not characteristic of a conventional Neel state. In- 
set shows the picture of single crystals used for diffraction ex- 
periments, (d) The temperature dependence of Bragg peak 
intensity of the (0, 1, 0.5) magnetic reflection from single crys- 
tal neutron diffraction measurement. The solid line is the 
power law fit described in the text. 

The systematic absences in the single crystal x-ray 
diffraction measurements unambiguously determine that 
the space group of Na 2 IrC>3 is C2/m and not C2/c as 
initially reported 



14] 



The typical crystal diffraction 
pattern shows diffuse streaking, characteristic of stack- 
ing faults within the layer sequence. Stacking faults in- 
volving fractional translation and rotation of the funda- 
mental C2/m layer module have been modeled to some 
extent in the isostructural Li2MnC>3 [l6j]. Polytypism 
analogous to that observed in the micas is also possible. 
We adopted a structural model that allowed for inter- 
mixing of the Nal and Ir sites to artificially account for 
some amount of stacking disorder, yet retain the ideal 
stoichiometry. The overall structure exhibits a virtually 
regular honeycomb layer of edge-sharing IrOg octahedra, 
similar to that observed in other so-called trioctahedral 
sheets [e.g., gibbsite Al(OH)3] in which the octahedra 



are slightly flattened perpendicular to layer stacking. In 
addition, the three O-Ir-0 bond angles perpendicular to 
the basal plane are all greater than 90° whereas the bond 
angles across the shared edges are 84.1(3)° and 84.5(3)°, 
in contrast to the undistorted 90°, as shown in Fig. 2(a). 
The structural distortion indicates a presence of the trig- 
onal crystal field in addition to the cubic crystal field, 
due to the repulsion of neighboring Ir atoms across the 
shared-edge of the octahedra. The trigonal crystal field 
in Na2lrC>3 makes the otherwise well separated gap be- 
tween J c ff' = 1/2 and 3/2 levels Q less discernible or 
pronounced, highlighting an important role in determin- 
ing the electronic band structure topology 17| . 



TABLE I. Structural parameters at T=125 K from single crys- 
tal x-ray diffraction measurements. The full data sets could 
be indexed using space group C2/m with a = 5.319(1) A, 
b = 9.215(2) A, c = 5.536(1) A and (3 = 108.67(1)°. The 
Ir-0 bond distances are 2.069(8), 2.067(9), and 2.060(12) A, 
and the Ir- ■ -Ir distances are 3.073(1) and 3.0705(8) A. Refine- 
ments are made using SHELXL-97 [Tq|. yielding an agreement 
factor Rl = 0.0687 for 334 reflections with F ob3 > 4a(F obs ). 
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The magnetic ground state is further characterized 
by the neutron diffraction on the single crystals. The 
magnetic propagation wave vector was determined to be 
q m = (0,1,0.5) in the C2/m notation based on exten- 
sive survey in reciprocal space using the four-circle neu- 
tron diffractometer. Figure 1(d) shows that the mag- 
netic Bragg peak intensity (Ib oc |M s | 2 , M s is the order 
parameter) disappears above Tn — 18.1 ± 0.2 K, con- 
sistent with the anomaly observed in the specific heat 
data (it is remarkable that the anomaly of C(T) down- 
shifts only by 0.5 K at 14 T.) [Fig. 1(c)]. Fitting I B to 
the power law scaling function of (1 — T/T)v) 2 ^ yields 
a critical exponent /3 = 0.29(2) that describes the mag- 
netic system to be three dimensional. The determina- 
tion of magnetic propagation wavevector and the correct 
description of the crystal structure put stringent con- 
straints on the possible magnetic models. The magnetic 
reflection appearing at (0,1,0.5) rules out the Neel config- 
uration [characterized by antifcrromagnetically-coupled 
nearest neighboring spins with qj n = (0,0,0.5)] but 
leaves the choice of either stripe or zigzag order in the 
basal plane as depicted in Figure 2(c)-2(d). Group the- 
ory analysis indicates that the magnetic representation 
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FIG. 2. (Color online) (a) Local structure within the basal 
plane. The compression of Ir06 octahedron (indicated by ar- 
rows) along the stacking leads to the decrease of O-Ir-0 bond 
angles across the shared edges, (b)-(c) Comparison of stripe 
and zigzag order that are consistent with the symmetry asso- 
ciated with observed magnetic reflections. In both cases, the 
Ir moments between honeycomb layers are antiferromagneti- 
cally coupled. 



r m ag can be decomposed into irreducible representation 
(IR) r mag = Ti + T 2 + 2r 3 + 2r 4 with corresponding 
basic vectors (BV) listed in Table II. Since the moment 
direction has been characterized to be along the a axis by 
magnetic susceptibility measurement and polarized x-ray 
15| . this information is implemented to per- 



technique 

form the model calculation and magnetic structural re- 
finement. Figs. 3(a)-3(c) show the rocking scans of three 
characteristic magnetic Bragg reflections qi = (0, 1,0.5), 
q2 = (0, 3, 1.5) and q3 = (0, 3, 0.5) in the (0, k, I) scatter- 
ing plane. The strongest reflection occurs at qi and the 
intensity decrease sharply at q 2 that has a larger momen- 
tum transfer. In contrast, there is no sign of magnetic 
scattering at q3 at base temperature. For single crystal 
magnetic scattering at wavevector transfer q, the mea- 
sured intensity follows: 



|^±(q)| 2 = |F m (q)| 2 -[e.F m (q)] : 



(1) 



where e is the unite vector along the q, and F m (q) is the 
magnetic structure factor that can be expressed as 



F m (q) = pV/j(q)Sk,, exp 27ri(q ■ rj). 



3=1 



(2) 



Here the sum is over all the magnetic atoms in the crys- 
tallographic cell, p = (r e 7/2) = 0.2695, Skj are the 
Fourier components proportional to the BVs listed in Ta- 
ble II, r is the vector position of atom ■?', a nd /(q) is the 
magnetic form factor for the Ir 4+ ions [19| . 

As summarized in Table III, both stripe and zigzag 
spin orders give the identical ratio |FL(q2)/Fi(qi)| 2 . 



TABLE II. Basis vectors (BVs) tpi of an IR of the space group 
C2/m and k = (0,1,0.5). BVs are defined relative to the 
crystalline axes. Magnetic moments for j atom in I th cell is 
given by m t j = ^ k S k j exp(-27rik-Ri) and S kj = dtpi, 
where d is the mixing coefficient. Only F3 and T4 are relevant 
since they describe the correct spin direction along the a-axis. 

^i(ri) ^ 2 (r 2 ) y> 3 ,^4(r 3 ) j> 5 ,M^i) 

Ir (0,0.333,0) (0,1,0) (0,1,0) (1,0,0), (0,0,1) (1,0,0), (0,0,1) 
Ir (0,0.667,0) (0,1,0) (0,-1,0) (1,0,0), (0,0,1) (-1,0,0), (0,0,-1) 



TABLE III. Calcuated magnetic scattering [Fj_(qi)[ 2 at qi 
(normalized to 100), q 2 and q 3 for stripe and zigzag spin 
orders and their comparison to the measurement. Errorbar is 
statistical and refers to one standard deviation. 





qi = (0,1,0.5) 


q 2 = (0,l,1.5) q 3 


= (0,3,0.5) 


T 3 (stripe) 


100 


51.1 


186.5 


T 4 (zigzag) 


100 


51.1 


0.002 


measurement 


5.40 ±0.38 


2.77 ±0.32 






Therefore the magnetic scattering at these two reflec- 
tions alone can not distinguish the difference between 
the two spin configurations. However, the magnetic scat- 
tering at q 3 is expected to be strong for the stripe spin 
configuration but absent for the zigzag spin chains; the 
absence of the magnetic scattering at q 3 illustrated in 
Fig. 3(c) clearly indicates a presence of the zigzag spin 
order. To determine the magnitude of the magnetic mo- 
ment/Ir, a large set of nuclear reflections under the same 
experimental configuration were collected to get the scale 
factor for normalization, yielding a magnetic moment 
of 0.22(1) /zb/R; this is considerably smaller than that 
(1 a*b/R) for a S — 1/2 system, consistent with early ob- 
servations for systems such as Sr2lr04 and Balr0 3 where 
the ordered moment is no more than 15% of 1 /x B /R 0- 
H, [2^ . The significantly reduced moment might be as- 
cribed to the strong hybridization of the Ir 5d orbital with 
the ligand oxygen 2p orbital and the moments are largely 
canceled out in the antiferromagnetic state. Moreover, 
the wavevector scans presented in Figs. 3(d)-3(e) show 
resolution limited Gaussian profile and Lorentzian type 
lineshape for in-plane and out-of-plane magnetic corre- 
lation, respectively. The data reinforce that the spins 
form a long-range ordered structure in the honeycomb 
basal plane while a short range order might still remain 
to some extent (with correlation length £ f» 139 ± 21 A) 
between layers due to the inherent imperfection in crys- 
tal structure as suggested in x-ray diffraction presented 
above. 

In recent theoretical proposals, Na2lr0 3 is regarded as 
one of a few model systems that can be mapped into 
the exactly solvable Kitaev model 2l|. The combina- 
tion of isotropic Heisenberg exchange interaction and 
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FIG. 3. (Color online) The rocking scans of characteristic 
magnetic reflections of (a) (0,1,0.5), (b) (0,1,1.5) and (c) 
(0,3,0.5). (d)-(e) The in-plane and out-of-plane wavevector 
scans for the (0, 1, 0.5) peak. Solid line in panel (e) is the fit 
to the Lorentzian form of A/(l + q 2 //i 2 ) (k being inverse cor- 
relation length) with instrument resolution convoluted. Hor- 
izontal bars denote the instrument resolution. 



anisotropic Kitaev term through strong spin-lattice cou- 
pling gives rise to a rich variety of the low energy mag- 
netic ground states. This includes the topologically non- 
trivial quantum spin Hall system in the weak interac- 
tion limit and evolution from the conventional Neel 
order to the spin liquid state sandwiched by a stripe 
phase depending on the microscopic parameters in the 



strong spin-orbit coupling limit [12j, |22|. The geomet- 
ric frustration due to the longer range exchange paths 
and the dynamic frustration caused by the Kitaev term 
leave the physical properties of Na2lrC>3 highly tunable 
by small perturbations such as magnetic field, vacan- 
cies and structural distortions 2^-26|. Only recently 
has the magnetic ground state been experimentally ex- 
amined and proposed to be a possible zigzag spin state 
using resonant magnetic x-ray scattering 15[ . The unex- 
pected spin state inconsistent with the original Kitaev- 
Heisenberg model underscores the novelty of the ground 
magnetic state, prompting theoretical suggestions that 
that the zigzag magnetic order could be explained only 
when the long-range magnetic Heisenberg interactions 




231 ] or a trigonal distortion of the IrOg octahedra 
in the [1,1,1] direction (local basis of the octa- 
hedron) is taken into account. Indeed, a new quantum 
phase transition from normal to topological insulator was 
recently predicted in Na2lrC>3 if both the long-range hop- 
ping and trigonal crystal field terms are included 17] . On 
the other hand, noticeable inconsistencies still exist in the 



band topology predictions that are likely due to the struc- 
tural parameters used for the first-principles calculation 



29j . With the presence of the trigonal crystal field, 
it is suggested that the J ff = 1/2 doublet is no longer 
as critical in Na2lr03 as in Sr 2 Ir04 dHl]) BalrOa @ and 
other iridates; instead, the trigonal crystal field (0.6 eV) 
and long-range hoping dictate the topological character, 
which is extremely sensitive to slight structural changes 
.17]. 

One of the unique aspects of this work is that both neu- 
tron scattering and x-ray diffraction data were collected 
from single crystals of Na2lr03. The results of this work 
therefore provide the much-needed, well-defined charac- 
teristic of the magnetic and crystal structures of the 
honeycomb lattice, thus significantly improve our under- 
standing of this intriguing system. We expect this study 
will clarify the topological character of the ground state 
in Na2lr03, a fertile ground yet to be fully explored. 

We thank S. Okamoto and C. de la Cruz for invaluable 
discussions. Research at the Oak Ridge National Labo- 
ratory was sponsored in part by the Division of Chemical 
Sciences, Geosciences, and Biosciences, and the Scientific 
User Facilities Division, Office of Basic Energy Sciences, 
U.S. Department of Energy. The work at University of 
Kentucky was supported by NSF through grants DMR- 
0856234 and EPS-0814194. 

Note added: after finishing our work, we become aware 
of recent inelastic neutron scattering work on the powder 
samples of Na 2 Ir03 and single crystal x-ray diffraction 
studies that reveal a similar structural information [30(. 
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